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We present atomistic calculations of quantum coherent electron transport through fulleropyrroli- 
dine terminated molecules bridging a graphene nanogap. We predict that three difficult problems in 
molecular electronics with single molecules may be solved by utilizing graphene contacts: (1) a back 
gate modulating the Fermi level in the graphene leads facilitate control of the device conductance in 
a transistor effect with high on/off current ratio; (2) the size mismatch between leads and molecule 
is avoided, in contrast to the traditional metal contacts; (3) as a consequence, distinct features in 
charge flow patterns throughout the device are directly detectable by scanning techniques. We show 
that moderate graphene edge disorder is unimportant for the transistor function. 

PACS numbers: 73.63.-b, 73.63.Rt, 85.65.+h, 72.80.Vp 



I. INTRODUCTION 

The isolation of graphen^ has lead to a wide range of 
scientific discoveries and technological opportunities.!^ 
These include, for instance, a half-integer quantum Hall 
effeclP^ with potential for a drastically improved quan- 
tum resistance standardj'^'and a high electron mobility in 
this atomically thin crystal for applications in radiofre- 
quency electronic^^ with reduced short channel effectPI 
or as a transparent flexible electrode.'^ In a wide per- 
spective, graphene is a material with potential as a versa- 
tile and controllable bridge between the atomic and the 
micron scales, with unique opportunities for nanoelec- 
tronics applications. Chemistry tools may alter graphene 
properites either globally (for example chemical ga tinjjS I) 
or locally (atoms and molecules binded to graphend^^E^ . 
At the same time, modern lithographic techniques com- 
patible with semi-conducting technology can be used to 
pattern graphene into devices and integrate them with 
conventional electronics .^^^ Building on these discoveries, 
we show in this paper how graphene patterned to form 
a nanogap can be used as electrodes for gate-tunable 
molecular electronics with single molecules. The transis- 
tor effect is achieved by gating the graphene electrodes 
while additional device functionalities can be built into 
the molecule bridging the nanogap. 

The idea of utilizing single molecules as active elements 
for electronics applications are based on a number of 
observations,1S"i^ including device minituarization, re- 
producibility, and functionality. Besides being the ulti- 
mately small object, molecules can be mass-replicated 
and their functionality can be tailored through molecu- 
lar synthesis. Many functional single molecule devices 
have been demonstrated,E^ but many problems remain 
before practical applications can be realized. Tradition- 
ally, a metal such as gold has been used to make con- 
tacts, although other configurations such as semicon- 
ducting substrates combined with the scanning tunneling 
microscope^*' (STM), have been shown to work as well. 
The huge size mismatch between metallic leads (> 10 
nm thick) and molecules is unavoidable, which makes 



nanogap fabrication challenging. Utilization of the STM 
for contacting molecules is not a scalable technology and 
mainly suited for making devices for research. In ad- 
dition to the problem of making nanogaps, an equally 
important problem for molecul ar ele ctronics is how to 
make good molecular transistors .^221221 xhe difficulty is to 
put a gate sufficiently close to the molecule in a metallic 
nanogap to achieve a gate effect. 

Here, we study theoretically the prospects of using 
graphene as a platform for single molecule electronics, 
where graphene nanostructures are used as contacts and 
interconnects instead of metal wires. The main purpose 
of this paper is to show that a transistor effect can be 
achieved by utilizing a back gate that changes the elec- 
tron density and the density of states at the Fermi level 
of the graphene leads. This transistor effect works well 
when the coupling of the gate to the molecule is weak 
compared with the coupling of the gate to the leads, 
which is the opposite situation to a traditional molecular 
transistor, including nanotube-based device^^H although 
sharing the same robustness as these through a resonance 
tunneling mechanism. In addition to the transistor ef- 
fect, the usage of graphene, being one-atom thick, would 
circumvent the size-mismatch problem experienced with 
metal contacts. 

The current fast improvements in graphene pattern- 
ing and device fabrication,^^ have opened new op- 
portunities for making advanced devices. The se inc lude 
gas sensors,^ nanopores for DNA sequencing,'22l3n! ^nd 
single-electron transistors operated as read out devices.fi^ 
In the latter experiment, magnetic molecules were de- 
posited on top of a graphene constriction. By utilizing an 
external magnetic field, the spin states of the molecules 
were manipulated, as could be read off by the graphene 
single-electron transistor working in its conducting state. 
We conclude that with this rapid progress in graphene 
device fabrication, the devices we shall study here can be 
made in the near future. 
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FIG. 1: (a) Geometry of the transistor, (b) _Diim6&e// molecule 
consisting of a 1,4-plienylenediamine bridge with C60 anchor- 
ing groups, (c) Energy level diagram at zero source-drain 
voltage. The work function W of the leads can be changed by 
a few hundred meV via a back gate so that the Dirac point 
Ed in the graphene density of states N{E) is above or below 
(depicted) the Fermi level Ef- The estimated charge transfer 
to the molecule, reflected in the contact potential CP, is such 
that the Fermi level intersects the LUMO. 



II. MODEL 

The geometry of the molecular transistor we are con- 
sidering is shown in Fig. [ija). The molecule, resem- 
bling a dumbbell, consists of a central wire of 1,4- 
phenylenediamine with C60 end-groups (i.e. fulleropy- 
rrolidine terminated benzene) , as depicted in Fig. [ijb) . 
The wide graphene leads extends far from the molecule 
and are electrically connected to source and drain. A 
back gate can be used to change the position of the Dirac 
point in the graphene bandstructure with respect to the 
molecular energy levels. We describe this gate effect in 
more detail below. 

The goal of this work is not to predict the functionality 
of a certain molecule in great detail, or to reproduce or 
explain a certain experiment. Rather, the goal is to show 
the most salient features of a single-molecule transistor 
with graphene leads operating in the quantum coherent 
regime. But we base our studies on the specific molecule 
in Fig. [T] This molecule has recently attracted a lot of 
attention because it has promising properties for mak- 
ing reproducible contacts via the large C60 anchoring 
groups. l^iiSlI Thus, both leads and anchoring groups are 
made of carbon. 



We shall in this work use a minimal model based on a 
tight-binding description. The Hamiltonian is 

where Ei are onsite energies, Uj are hopping amplitudes 
between sites j and i, and the operators c] and Ci create 
and destroy electrons on sites i. We concentrate on the 
quantum coherent transport regime and leave effects of 
Coulomb interaction to future studies. This corresponds 
to assuming a small charging energy U on the molecule 
compared with the molecular level broadening F due to 
the coupling to the leads. The graphene nanostructured 
leads, as well as the molecule, are readily built up by re- 
stricting the sites i and j. We study both armchair and 
zigzag graphene nanoribbon leads with nearest neight- 
bor hopping amplitude t, with and without edge disor- 
der. We include large and wide sections of the graphene 
leads in the calculations and connect them to ideal rib- 
bons connected to reservoirs through the technique of 
self-energies. This is a Landa uer ap proach based on non- 
equlibrium Green's functions .E^EH foj- simplicity 

focus on the low-bias regime and present results for the 
transmission function of the device, as well as spectral 
charge current flow patterns inside the device. We note 
that it is important to have wide ribbons, with edge dis- 
order on a length scale smaller than the ribbon width, 
otherwise weak links may form at necks of the imperfect 
graphene ribbon. Two such necks define a quantum dot 
with single-electron transistor properties,!^ which leads 
to unwanted Coulomb blockade effects in the leads. 

For the molecule, the C60 end groups and the ben- 
zene ring of the bridge are all carbon based, while each 
link group in addition contains a nitrogen atom. For our 
purposes, it is enough to model the molecule within the 
tight-binding (Hiickel) theory on equal footing with the 
leads, and leave details to be explored in future calcula- 
tions and experiments. The parameters of the molecule 
are kept to a minimum by only varying hopping between 
C60 end groups and the center phenylenediamine bridge, 
while letting all other carbon atoms have the same pa- 
rameters in the molecule as in the graphene leads. The 
molecule parameters are the onsite energy on the ni- 
trogen atom relative to the onsite energy of the carbon 
atoms,l^i?Ar = Ec- 0.9|t|, {t is the C-C hopping ampli- 
tude), the hopping from C60 to the nitrogen group 0.4i, 
and the hopping from the benzene ring to the nitrogen 
group 0.6t, see Fig.[2jb). The energy levels and orbitals of 
the C60 end groups are known within the tight-binding 
model.'^fil For the bridge, we have compared our model 
with orbitals obtained with the freely available quantum 
chemistry package GAMESS,'22l see Fig. [2];a). We would 
like to emphasize that if the parameters of this model are 
varied, unimportant details of the results may change, 
but the general principles of how the transistor operates 
will not change. 

Finally, to determine the Cgo-on-graphene adsorption 
geometry and estimate both the magnitude of effective 
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FIG. 2: (a) Molecular eigenvalues of the phenylenediamine 
bridge without C60 end-groups (black squares), the C60 
molecule (red rings), and the dumbbell molecule (blue stars) 
within the tight-binding model with the indicated parame- 
ters. The dumbbell HOMO originates from the phenylene- 
diamine bridge, while the LUMO originates from the end- 
groups. The insets show the HOMO (left) and LUMO (right) 
of the phenylenediamine bridge computed with the quantum 
chemistry package GAMESS.'^'^ The symmetries of these or- 
bitals are also obtained within the tight-binding model, (b) 
Geometry of the molecular bridge. Carbon atoms number 
1, 2, 7, and 8 are in sp"^ orbital hybridization and do not 
participate in current transport, but we include next-nearest 
neighbor coupling tcN = OAt between atoms 3, 4, 5, and 6 to 
the nitrogen atom (number 9). The hopping from nitrogen to 
the nearest neighbor in the benzene ring (atom number 10) 
is tjvc = 0-6t. (c) Bernal stacking of a hexagon in C60 (red) 
and graphene (black); orientation 02 in Fig. |3] (d) Stacking 
for a 30° rotation of the C60; orientation Ol in Fig. [s] The 
nearest neighbor hopping [vertical in (c) and dotted in (d)] 
and next-nearest neighbor hopping (dashed lines) between the 
two layers are ti = 0.28t and t2 = 0.22t, respectively. 



Cgo-graphene hopping constants and the Cgo-graphene 
charge transfer, we have performed a study usin g the 
van der Waals density functional (vdW-DF) method.'^SEni 
The results have been obtained in a non-selfconsistent 
evaluation of the most recent version, vdW-DF2,'^a ver- 
sion which some of us have previously found gives an 
accurate description of the binding in both a Cgo crystal 
and of graphene layers .132 We have used a plane- wave code 
and a standard semi-local density functional approxima- 
tion, to obtain underlying results for the electron-density 
variation (as a function of adsorption geometries and dis- 
tances). Our choice of non-selfconsistent vdW-DF eval- 



uations is motivated by a recent analysis.!^ 

We present in Fig. [3] a summary of our vdW-DF study. 
Panel (a) shows vdW-DF2 results for the variation in 
binding energy of Cgo on graphene with center-of-mass 
separation for three typical low-energy adsorption geome- 
tries that we have investigated in an extended search. We 
find that there is a systematic preference for adsorption 
with the Geo hexagon facing down and situated on top of 
a graphene atom. The panel also provides a comparison 
of this type of adsorption geometries (identified by the in- 
serts which shows Cgo atoms in purple, graphene atoms 
in black) and we find that the most favorable configura- 
tion corresponds to a 30° rotation of what would amont 
to a Bernal stacking of the hexagon on the graphene. Im- 
portant for the transport modeling, we predict that the 
optimal adsorption separation is smaller than the value 
(vertical dashed line) which would corresponds to the 
predicted layer separation in graphite. We conclude that 
the effective Cgo-to-graphene hopping constants must be 
chosen larger than the choice which is made in a tight- 
binding modeling of graphite; in the qualitative transport 
modeling below, we simply set the enhancement at a fac- 
tor of two. The two lowest energy configurations (solid 
and dashed lines in Fig. |3]) corresponds, in the transport 
calculation below, to anchoring of the dumbbell molecule 
for zigzag (orientation 01) and armchair (orientation 02) 
leads, respectively [see Fig. [2|^c)-(d)]. 

In Fig. [Sj^b) we show the details of the Ceo adsorption 
and the complex charge transfer which we have calcu- 
lated for the optimal adsorption geometry. The blue col- 
ors shows regions of electron accumulation whereas the 
red regions identify an electron depletion (relative to a su- 
perposition of the graphene and Cgo electron densities). 
The panel shows that the binding is characterized not 
only by van der Waals forces but also by a pronounced 
dipole (and even multipole) formation. In addition, we 
fincP^ a net charge transfer from the graphene and to the 
60- We find that the Cgo-on-graphene binding is be- 
yond simple physisorption as the charge rearrangement 
causes a work function shift .'^ESI 

The first step in the transport characterization is to 
obtain the retarded Green's function G^{E) of the sys- 
tem, which is a matrix in the site indices. We utilize 
our own implementation of a recently developed recur- 
sive algorithnPS within which the sites are added one by 
one, which is ideal for our devices with complicated ge- 
ometries. The advanced Green's function is obtained by 
hermitian conjug ation C^iE) = [G"iE)]\ Obscrvables 
are related to the lesser Green's function G^. In the 
absence of electron correlations, the expression for the 
lesser Green's function is reduced to the form 

e 

xJ2GgiE){[j:fiE)]^ - J:fiE)]^^Gi^{E). 

cc 

It involves the distribution functions of the leads fi{E) 
and self-energies I]f {E) at the surfaces of the leads that 
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remain after eliminating the leads in favor of the system 
Green's function. The leads are enumerated by the index 
I (here l — \ and 2 for source and drain) and surface sites 
of the leads are labeled by c and c. Local charge current 
flow in the device (bond current between sites i and j) is 
written as 

/oo 
{U,G<{E)~t,,G<{E)\dE. (2) 
-oo 

The transmission function can also be written in terms 
of the retarded Green's function and self-energies of the 
leads, 

nE) ^ Tr ^,{E)G'^^{E)Y^{E)Gi^{E)\ (3) 

where Vi = z[S^ — (S^)^'], Gfj is the propagator between 
leads 1 and 2, and the trace is over the surface sites. Or 
we may compute T{E) by integrating the bond-currents 
[Eq. [2] flowing through an interface of the device. 




FIG. 3: Results of a density functional theory calculation in- 
cluding the van der Waals interaction, (a) Binding energy 
of C60 on graphene for three high-symmetry configurations. 
The solid black line (orientation Ol) and dashed black line 
(orientation 02) correspond to dumbbell molecule anchoring 
for zigzag and armchair leads, respectively. In both cases, a 
hexagon of the C60 faces graphene. The configuration with 
a pentagon facing graphene (red line) is less favorable, (b) 
The charge density distribution at the C60-graphene contact 
for orientation 01. The blue (red) area is negative (posi- 
tive) charge, implying a visible charge transfer to C60 from 
graphene. 



III. TRANSISTOR EFFECT 

In Fig. |4|^a) we show an example of a transmission 
function for one molecule in the center of the graphene 
nanogap in a symmetric position, here for armchair leads. 
The transmission displays typical resonance features near 
the molecular levels of the isolated molecule. The levels 
are shifted and broadened by the coupling to the leads. 
The amount of broadening depends on the exact coupling 
of the molecule to the graphene, but also on the nature of 
the molecular orbitals. For this particular molecule, the 
LUMO is mainly centered on the C60 anchoring groups 
that act as effective extensions of the leads, while the 
bridge acts as the weak link in the system. Functional- 
ity can be added to the device by choosing a different 
bridge by exchanging the benzene ring during molecular 
synthesis."^ But here we shall continue working with the 
benzene bridge and focus on the transistor effect. 

In Fig. |4|^b) we present a contour plot of the transmis- 
sion function for energies (vertical axis) near the LUMO 
as we rigidly shift the band structure of the graphene 
leads relative the molecular level by a back gate volt- 
age (horizontal axis with a transfer function a between 
the gate voltage and the shift of the Dirac point). The 
gate effect we have in mind can be visualized start- 
ing from Fig. [T] as moving the graphene bands verti- 
cally keeping molecular levels and the Fermi level fixed. 
When the Dirac point is far from the molecular level, 
the transmission resonance corresponds to the transistor 
in the on-state. As the Dirac point in the bandstruc- 
ture approaches the molecular level, either from below 
or above in energy, the transmission resonance is shifted 
due to hybridization with zigzag nanogap edge states. 
The possibility of such hybridization was also noted re- 
cently in a DFT calculation.''^. When the Dirac point 
passes through the level, the transmission is suppressed 
and the transistor is in the off-state, see Fig. Qb) at 
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aeVg = i?LUMO ~ 0.133t. The on and off states will be 
well separated when the Dirac point can be shifted by 
dE]j > r, larger than the molecular level broadening F. 

The transistor on-off ratio will be large when the Fermi 
level Ep is aligned close to a molecular level. This is the 
case here, with Ep in the broadened LUMO. We can 
estimate the molecular level alignment^ with respect to 
the Fermi level by estimating the charge transfer between 
graphene and C60. We do that by comparing the work 
function of graphene with metals for which the charge 
transfer to C60 has been measured. It has been shown 
by scanning Kelvin probe microscopjH^ that the appli- 
cation of a back gate voltage results in a change of the 
graphene work function between 4.5 eV (electron doped) 
and 4.8 eV (hole doped). Scanning tunneling experi- 
ments and DFT calculations of C60 on gold and silver 
surfaces sho'w^Sl that the charge transfer to C60 from 
gold is vanishingly small, while it is of order 0.2e from 
silver. The work function of silver is 4.6 eV, while that 
of gold is 5.3 eV. This picture is corroborated by our 
vdW-DF study which identified a net charge transfer and 
dipole formation. In summary, we draw the conclusion 
that there is considerable charge-transfer effects already 
for hole-doped leads {W — 4.8 eV), which results in the 
Fermi level energy in our system aligned inside the broad- 
ened LUMO. As the gate voltage is changed, the charge 
transfer to the molecule will increase as we go through 
the Dirac (neutrality) point of graphene to the electron 
doped side (eventually reaching W = 4.5 eV), which re- 
sults in a Fermi level deeper inside the LUMO. Based on 
these estimates, where the Dirac point can be shifted by 
at least O.lt through the LUMO, see Fig.|4];d), the on-off 
ratio will be large but the precise value will ultimately 
have to be determined by experiment. The back gate is 
straightforward to use compared with the traditional di- 
rect gating of the molecule itself. In fact, the transistor 
effect is most effective when the gate is decoupled from 
the molecule and only affects the graphene leads. 

In Fig.[5ja) we show the gate effect for the case of wide 
zigzag graphene leads, with the anchoring of the C60 end 
groups in orientation 01. For this orientation, there is 
no nanogap edge states, since the nanogap has armchair 
orientation. The strong hybridization of the molecular 
level with lead states is therefor absent and we predict a 
simple weak shift of the molecular level with gate voltage. 
As the Dirac point is tuned through the molecular level, 
the transmission is quenched, and we have a transistor 
effect analogous to the case with armchair leads discussed 
above. 

We note that graphene itself (without nanogap and 
molecules as weak links) works as a transistor via the 
back gate. However, the graphene transistor can not be 
set in the off-state as the minimal conductivity is of or- 
der e^/h, in contrast to the molecular transistor with 
graphene leads that we study here. A nanoribbon has 
an energy gap related to its width, and would work as a 
transistor with an off-state. However, the required rib- 
bon width is small (a few nm) and it is very difficult 
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FIG. 4: Electron transport through a single dumbbell 
molecule connected to armchair graphene leads with perfect 
edges. Anchoring of the dumbbell molecule in orientation 02. 
(a) Transmission as function of energy when the Dirac point 
(here S = 0) is below the LUMO, -EluMO ~ 0.133i. (b) 
Transmission as function of energy and back-gate voltage for 
energies close to the LUMO. As the Dirac point is tuned by 
the back-gate through the molecular level, the transmission 
is quenched which leads to a transistor effect, (c) A sketch of 
the movement of the Dirac point through the molecular level 
as the graphene lead work function is tuned by the back gate 
voltage. We estimate (see text) that the charge transfer to 
the molecule increases as the Dirac point Ed in the graphene 
band structure is tuned to be below the Fermi level. 
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FIG. 5: Electron transport through a single dumbbell 
molecule connected to graphene leads, (a) Zigzag leads with 
perfect edges; anchoring of the dumbbell molecule in orien- 
tation Ol. (b) Armchair leads with random rough edges; 
compare Fig. Bib) for perfect edges. 
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to control the nanopatterning with the required atomic 
resolution. In contrast, as we show below, the molecular 
graphene nano-gap transistor is more robust against edge 
disorder in the graphene leads. 

IV. CHARGE FLOW PATTERNS 

Since graphene is 100 % surface, it is an ideal ma- 
terial to study by scanning techniques.'^ Previously, 
atoms and molecules on metal surface^^ or large-area 
graphen^Sfi! have been studied by scanning tunneling mi- 
croscopy (STM) and spectroscopy (STS) and valuable 
details about e.g. orbitals and charge transfer have been 
obtained. Scanning techniques have been utilized to 
reveal local information about Coulomb interactions in 
graphene nanostructures.l^ Quantum transport through 
quantum point contacts in 2DEGs have been mapped 
oulp2l by scanning gate spectrocsopy (SGS) and revealed 
so-called branched flow originating from a background 
random potential due to doping impurities in the lay- 



FIG. 6: (a) Spectral current flow pattern through a transistor 
with one molecule in the nanogap. Note that the molecule is 
much smaller than the lobe structured pattern in the charge 
flow in and out of the molecule visible in the leads, (b) Spec- 
tral current flow pattern for the case of two molecules in the 
nanogap. (c) Spectral current flow patterns disturbed by ran- 
dom edges of the leads. In all cases, we assume zero temper- 
ature and study the linear low-voltage regime. 

ers forming the 2DEG. SGS has also been used to re- 
veal cohere nt electron transport in large-area graphene 
flakes.'S^^ Clearly, transport in a molecular device with 
metal electrodes is hidden by the bulky nature of the 
metallic contacts. Graphene on the other hand, be- 
ing two-dimensional, would be a uniquely suitable elec- 
trode enabling information about quantum transport in 
a molecular device to be revealed by scanning techniques. 

In Fig. we present the spectral current flow pat- 
tern [the integrand in Eq. w through the device for an 
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energy corresponding to the top of the LUMO peak in 
the transmission shown in Fig. |4][a) . The position of the 
molecule is clearly visible in the current flow pattern. 
The molecule forms the weak link where all current is 
channeled through. The current is flowing in and out of 
the molecule in a characteristic lobe pattern, that is due 
to the specific anchoring of the C60 on graphene. Deep 
inside the electrodes, the current is carried throughout 
the width of the ribbon. 

In Fig. ^et) the molecule is in a symmetric position in 
the nanogap. If the molecule is in an asymmetric posi- 
tion, the charge flow pattern is simply displaced verti- 
cally and changes in an intuitive way (not shown) and 
the transmission function remains unchanged unless the 
molecule is very close to the upper or lower edges of 
the nanogap, within a few rings in the graphene leads. 
Similarly, if two molecules are bridging the nanogap, to 
a good approximation two flow patterns are simply su- 
perimposed, as we show in Fig. |6jb). We may expect 
quantum intereference between the two pathways (two 
molecules). The effect on the transmission function is 
weak in this example, however, unless the two molecules 
are very close to each other, with the anchoring groups 
only separated by a few rings in the graphene leads. 

In a real device, perfect armchair or zigzag edges are 
hard to fabricate. In Fig. ^c) we show an example of the 
effect of edge disorder, consisting of randomly removed 
carbon atoms within one ring from the original perfect 
edges. The interference patterns in the leads are now 



affected, but the molecular weak link remains clearly vis- 
ible in the patterns of current flowing in to and out of 
the molecule. Also the transmission function T{E) is in 
its main features unaffected by defects in the leads, as 
we show in Fig.jsf^b) [compare Fig.|4|^b)]. Ideal graphene 
leads are not crucial for the transistor to function, since 
the weak link is the molecule. 



V. SUMMARY AND CONCLUSIONS 

In conclusion, we have studied a single molecule device 
with graphene leads working in the quantum coherent 
transport regime. We predict a transistor effect that is 
pronounced when the gate coupling to the leads is much 
stronger than to the molecule. This opens new avenues 
for research of gate tunable quantum coherent molecular 
electronics with single molecules. 
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